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^ , In this chapter we propose an extended concept of the expected discounted penalty function (EDPF) 

' that takes into account new ruin-related random variables. We add to the EDPF, which was introduced in 

classical papers [Gerber and Shiu (1997), (1998) and Gerber and Landry (1998)|, a sequence of expected dis- 
counted functions of new record minima reached by a jump of the risk process after ruin. Inspired by results 
of Huzak et al. (2004) and developpements in fluctuation theory for spectrally negative Levy processes, we 
provide a characterization for this extended EDPF in a setting involving a cumulative claims modelled by a 
subordinator, and Brownian perturbation. We illustrate how the extended EDPF can be used to compute 
' the expected discounted value of capital injections (EDVCI) for Brownian perturbed risk model. 

Keywords: ruin, spectrally negative Levy process, scale function, Gerber-Shiu function, Laplace transform, 
capital injections. 
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^ ! 1 Introduction 

lO ■ The concept of Expected Discounted Penalty Function (EDPF) has been introduced in classical 

! papers [Gerber and Shiu (1997), (1998)]. This so-called Gerber-Shiu function is a functional 

of the ruin time (i.e., the first time the reserve level of a firm becomes negative), the surplus 
prior to ruin, and the deficit at ruin. The EDPF has been extensively studied and generalized 
to various scenarios and there is now a wide range of models for which expressions of the EDPF 
are available. Since the EDPF operates on a random cashflow at ruin, where the cashflow is 
^ ■ a function of the deficit at ruin and the surplus prior to ruin, applications in the context of 

, insurance and finance are quite natural. For example, the EDPF can be used to determine the 

initial capital required by an insurance company to avoid insolvency with a minimum level of 
confidence and for fixed penalization of the ruin event. Similary, the EDPF can be used as a 
pricing device for American options [Gerber and Shiu (1998b)]. 

In oder to have a more valuable Gerber-Shiu function for the management of insurance 
risks and the monitoring of the solvency of a firm, Biffis and Morales (2010) extended the 
EDPF to include path-dependent penalties. In particular, they generalized the definition of 
EDPF to include a new random variable, the last minimum of the surplus before ruin. They 
obtained a defective renewal equation for this generalized EDPF. The representation is obtained 
for a subordinator risk model perturbed by a spectrally negative Levy process. More generally, 
when the risk process is driven by a spectrally negative Levy process, Biffis and Kyprianou 
(2010) provided an explicit characterization of this generalized EDPF in terms of scale functions. 
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extending results available in literature. One of the reasons for the limited use of such EDPF 
is that the last minimum of the surplus before ruin, the surplus prior to ruin and the deficit 
at ruin, only characterize the surplus before and in a neighborhood of the ruin time. In other 
words, none of the arguments in the EDPF can be used as a predictive tool for successive deficit 
times after ruin. The situation would change if a penalty could apply after ruin, for example 
by acting on relevant characteristics of the paths of the risk process that may lead to successive 
minima after ruin, and not just on its level before, at, or immediately prior to ruin. 

In this chapter we show how to extend the EDPF to include these new random variables. In 
particular, we generalize the EDPF to include the sequence of successive record minima reached 
by a jump of the risk process after ruin. We obtain a new form of EDPF which gives character- 
istics of the paths of the risk process after ruin and not only before, and in a neighborhood of, 
ruin time. There are practical applications of this extended EDPF in the context of insurance 
and reinsurance. For example, it can be used to determine the capital required by an insurance 
company to survive not only in the neighborhood of ruin, but also after ruin when the risk 
process continuous to jump downwards, that is it continues to pay out claims. Similary, this 
extended EDPF can be used by insurance company or government institutions to determine 
the capital which should be injected at each deficit time at and after ruin that will allow it to 
continue its operations. At the time of ruin, the insurance company could have access to other 
reserves allowing it to survive and pay to the customers claims made after ruin. It is at this 
time that the company might need the capital injections that will allow the net aggregate cash 
inflow to return to pre-ruin levels. This ruin cycle may occur several times and requires repeated 
interventions with capital injections. Thus, it would be interesting to add to the EDPF the ex- 
pectation of a sequence of discounted functions of the successive minima reached by claims of the 
risk process after ruin. This extended EDPF is needed to determine the Expected Discounted 
Value of Capital Injection (EDVCI) for a subordinator risk model perturbed by a Brownian mo- 
tion. Moreover, our approach provides a connection with some of the works in Einsenberg and 
Schmidli (2011) where they have studied the EDVCI for the classical risk model. Consequently, 
we generalize some results in Einsenberg and Schmidli (2011) where a similar problem is solved 
for the classical risk model. 

We use the results of Huzak et al. (2004) and developments in fluctuation theory for spec- 
trally negative Levy processes to give an explicit characterization of this extended EDPF. The 
characterization is obtained for a subordinator risk model perturbed by a Brownian motion. 

The chapter is organized as follows. In Section 3.2, we describe a perturbed subordinator risk 
model and define our extended EDPF. In Section 3.3 and we review and provide some results that 
are needed in our derivations. In particular, we briefly review some results in Huzak et al. (2004) 
about fluctuation theory for spectrally negative Levy processes, and we give some preliminary 
results for first-passage times of a subordinator risk process perturbed by a Brownian motion 
under a change of measure. In section 4, we provide the expression of the extended EDPF in 
terms of convolution product of densities which are identified in Section 2. In section 5, we show 
how the results of section 4 can be used to give explicitly the EDVCI for a subordinator risk 
model perturbed by a Brownian motion, and how it can be used to recuperate the expression of 
the EDVCI for the classical risk model [see Einsenberg and Schmidli (2011)]. Finally, Section 6 
offers some concluding remarks. 

2 Risk model and the Expected Discounted Penalty Function 

Let (ri,F,P) be a filtered probability space on which all random variables will be defined. Let 
us define S = {St)t>o to be a subordinator (i.e., a Levy process of bounded variation and 
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nondecreasing paths) without a drift. Let v be the Levy measure of 5; that is, is a cr-finite 
measure on (0, oo) satisfying ^^^^-^^^ l\y)v(dy) < oo. 

We define the spectrahy negative (i.e. a Levy process with negative jumps) process X = 
{Xt)t>o as 

Xt = ct-St + Zt, (2.1) 
where Z is a multiple of a standard Brownian motion, we write 

Zt = aBt, 

with B a standard Brownian motion independent of S. Let T := {J~t)t>o be the filtration 
obtained by a{Ss, Bs,s > 0). 

We consider a very general setup that generalizes the standard Cramer-Lundberg model. 
The model discussed in this paper is, 

Rt:=x-Yt, t>0, (2.2) 

where Y = (lt)t>o is a spectrally positive Levy process (i.e. a Levy process with positive jumps) 
defined by 1^ = —Xt, where Xf is given by (12. ip . As introduced previously, note that the risk 
process given by 

Rt = x + ct-St + Zt, (2.3) 

is on the same spirit as the original perturbed model in Dufresne and Gerber (1991). The 
constant x > represents the initial surplus, while the process Y represents the cash outfiow 
of an insurance company. The subordinator S represents cumulated claims, and this is why 
we need it to be increasing since the jumps represent claims paid out. The Brownian motion 
Z accounts for any fiuctuations affecting the components of the risk process dynamics, such as 
claims arrivals, premium income and investment returns; ct represents premium inflow over the 
interval of time [0,t]. 

The premium rate c is assumed to satisfy the net profit condition, precisely £[6*1] < c, which 
requires 

[ yv{dy) < c . (2.4) 

J(0,oo) 

The condition in equation (12. 4p implies that the process Y has a negative drift in order to avoid 
the possibility that R becomes negative almost surely. This condition is often expressed in terms 
of a safety loading. Indeed, it is standard to write the drift component within y as a loaded 
premium. For instance, notice that we can recuperate the classical Cramer-Lundberg model if 
£7 = where c := (1 + 0)E[S'i] and S is a compound Poisson process modeling aggregate claims. 
The drift c, with a positive safety loading > 0, is the collected premium rate. We do not use 
the concept of safety loading in this paper in order to simplify the notation but we stress the fact 
that this concept is implicitly considered within the drift of Y when we impose the condition in 
equation (|2.4p . The classical compound Poisson model can be incorporated in this framework by 
setting v{dy) = XK{dy), where A is the Poisson arrival rate and is a diffuse claim distribution. 

We refer to Asmussen (2000) for an account on the classical risk model, and to Furrer and 
Schmidh (1994), Yang and Zhang (2001), Huzak et al. (2004) and Biffis and Morales (2010) for 
different generalizations and studies of model (12. 2p . 

Now, one of the main objects of interest in ruin theory is the ruin time, r^, representing the 
first passage time of Rt below zero when Rq = x, i.e. 

Tx := inf{t > : Yt > x} , (2.5) 
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where we set = +00 if i?^ > for all t > 0. 

Associated with the ruin time r^;, we have at least two other quantities that contain relevant 
information on the ruin event from a risk management perspective, namely the deficit at ruin 
—Rtx = — X and the surplus immediately prior ruin Rt^- = x — y_- . 

Gerber and Shiu (1998) studied the ruin event in the compound Poisson case by analyzing 
the joint law of all these quantities in one single object, the EDPF. In the following, we define 
the EDPF under the model (lOl. 



Definition 2.1 Let w be a non-negative Borel-measurahle function on xM+ such that w{., 0) 
0. For q >0, the EDPF associated with the risk process ^2. 2\) is defined as 



(j){w;x;q) =K e-'^'^^w{x -Y^-,Yr^ - x) . (2.6) 

Note that the condition w{-,0) = excludes the event {Yr^ = x}. This possibility is known 
as creeping and we chose not to consider it in our analysis. For simplicity, we assume that the 
function w assigns a zero penalty when ruin occurs by continuously crossing over zero. Notice 
that for a model like (j2.2p . this only happens when ruin is caused by the Brownian motion 
component of the process Y. 

Following the same order of ideas, we study the EDPF under the general context which 
gives relevant informations on and after the ruin event. More precisely, we generalize the EDPF 
defined in (j2.6p to include the quantities associated with the ruin time Tx and times sequence of 
successive minima reached by a claim of the risk process (12. 2p after ruin. This implies that some 
of notations related to record minima need to be introduced . 

Thus, we must define the first new record time of the running supremum 

T:=mi{t>0,Yt>Yt-} , (2.7) 

and the sequence of times corresponding to new records of Y reached by a jump of S after the 
ruin time. More precisely, let 



r(^) := Tx, (2.8) 



and inductively on {r^") < 00} 



r("+i) := m 



f{t>T^^\Yt>Yt-}. (2.9) 



Recall from Theorem 4.1 of Huzak et al. (2004) that the sequence (t(")).„>i is discrete, and, 
in particular, neither time non any other time is an accumulation point of those times. More 
precisely, r > a.s. and r*-"") < r^"^"^^ a.s. on {r^") < 00}. As a consequence, we can order the 
sequence (t^"'^)„>i of times when a new supremum is reached by a jump of a subordinator as 
< t(i) < r(2) <" . . . a.s. 

Let us introduce the random number N given by 

N := max{n : r^") < 00}, (2.10) 

which represents the number of new records reached by a claim of the risk process (|2.2p . In the 
following, we study the EDPF in a new context involving the deficits at times (r^"'^)„>i. More 
precisely, in this paper we set out to study the following extended EDPF for the model (|2.2p 
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Definition 2.2 Let F = {Fn)n>o be a sequence of non-negative measurable functions from R-|_ X 
]R+ to M, X and q > 0. The discounted penalty associated with the risk process \2.2i) . F and q is 
defined as 

N 

=E[5^e-'?^^"'F„(y,(„_i),y,(„));T, < cx)] . (2.11) 

n=l 

We assume in the previous definition that t^'^^ = r~ and 

Fi(-,x)=0. (2.12) 

Note that the condition given by (12. 12^ is used to exclude from calculation the event {1^^ = x}. 
P{F,q,x) is an extension of the classical EDPF defined in (|2.6p . In particular, it reduces to 
4){w; x; q) since we suppose that Fi{u, v) = w{x — u,v — x) and F„ = for n > 2. 



3 Preliminary results 

In this section, we will give some preliminary results for first-passage times of the risk process 
defined in (12.2^ under change of measure. These results are based on the works in Huzak et al. 
(2004) where the ruin probability has been studied for a subordinator risk model perturbed by a 
spectrally negative Levy process. This allows us to give a more detailed analysis of the extended 
EDPF defined in (I2TTD . 

Recall from Section [2] that S" is a subordinator defined on the filtred probability space 
(r2,F,P). The Laplace exponent of S is defined by 



M») = I [e"^ - i]y{dy) , (3.1) 

J(0,oo) 



where 
Note that 



E[exp(QSt)] = exp(tV5(a)) ■ (3.2) 
E[Si] = ^Jj(0+) = / yu{dy) = [ Hy,(x^)dy , (3.3) 

J(0,oo) JO 

where the last equality follows from integration by parts. As explained before, we assume 
throughout that M[Si] < oo. 

The Laplace exponent ^ of X defined in (j2.ip is defined by the relation 



E[exp(/3Xt)] = exp(tVx(/3)), 

(3.4) 



where 



^Px{P) = c/3 + ^s(-^) + Vz(/3) 

= ^5(-/3) + Z(/3) /3>0, 

where Zt = ct + Z and ip^{/3) = c/3 + ipz{P)- The last equality is due to the independence of 
S and Z. We refer to Bertoin (1996), Sato (1999) and Kyprianou (2006) for a comprehensive 
account on Levy process. 
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Let us introduce the distribution function G of — inft>o(^t) = supj>q(— ct — Zt). Using a 
method similar to Yang and Zhang (2001) (see also Huzak et al. (2004)), the Laplace transform 
of G can be shown to be given by 



oo 



G(/3) = / e-^yG{dy) 



c/3 



V'z(/3) 



(3.5) 



That is, 



G(/3) = — (3.6) 
C/3 + 



2 

since Zt = aBt is a multiple of a standard Brownian motion. Then, G is given explicitly as an 
exponential distribution function with parameter 2c/cr^, i.e. G has density 

G{dy) = ^e-^^dy . (3.7) 

We also introduce the parameter 



£72 



P := ^ = 1/ yuidy)eiO,l). (3.. 



c 



(0,oo) 



We denote by Y the supremum given by Yt = sup^x Yg, for t > 0. 

Let us introduce the density of the overshoot at time r which is defined by 

H{du) = F{Yr - Y^- e du; T <oo) , (3.9) 

where u > 0. We shall give in the following proposition the density H{-) in terms of Levy 
measure and the premium rate. 

Proposition 3.1 Let Y the spectrally-positive Levy process defined in \2.3\) . 

1- The distribution ofYj- — Y^- on the set t < oo is given by 

1 

H{du) = - u{du + y)dy; u>0. (3.10) 
c Jo 

2- The distribution ofYj- is given by 

F{Yr e du; T < oo) = H * G{du) ; n>0, (3.11) 
where H * G(-) denotes the convolution of H{ ) with G{-) defined by 



f{u)H * G{du) = f{y + v)H{dy)G{dv) , 

for all Borel set ^ o/M x M. 

Proof. 1- Let us suppose that / is a nonegative bounded Borel function. We prove firstly that 

E[/(y^ - F,-); r < cx)] = E[ T J{Yt - Yt)dt\ , (3.12) 

Jo 
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where /(y) = /(g f{u — y)^{u>y}'^{du). The proof of (13.120 fohows by an application of the 
compensation formula [see Bertoin(1996) p. 9 or Theorem 4.4. of Kyprianou (2006)] appUed to 
the Poisson random measure, with intensity measure dt , v{du) associated to the jump of Y . We 
have 



E[/(y.-y,-); T<oo] = nY.f^^t-Yt-)\y^^y^_^^^^^] 

t>0 
poo f 

Jo J\0.oo) ' ~ 



El 



J[0,oo) 

f{u-{Yt-Yt))l^^^y^_y^]dtu{du) 



'0 J(0,oo) 

= E[ r J{Yt-Yt)dt]. (3.13) 
Jo 

Hence, from Proposition 4.3. of Huzak et al. (2004), the expected occupation time measure of 
(I3.13P is given by 

^[[m-YMt] = ^±^Jjmy. (3.14) 

Since P(r = oo) = 1 - P(t < oo) = 1 - p [see Corohary 4.5. of Huzak et al. (2004)], ([31^ is 
equal to 

f{u)du = - / f{u - y)l{u>y}v{du)dy 

" J(0,oo) 



c - E[Si] 70 c 

1 
c 



oo 

f{u) I iy{du + y)dy. (3.15) 

(O.oo) 



Equating the left-hand side of (|3.12p and (|3.15|) implies that 

F{Yr - Y^- edu; T <oo) = - [ [ v{du + y)dy , (3.16) 

C J(0,oo) Jo 

and the statement 1 of Proposition 13.11 follows . 

2- Using the conditional independence of 1^ — Y^- and Y^- given r < cxd, 

P(yr e dti ; T < oo) = P(y^ G dnj r < oo)P(r < oo) (3.17) 
= P(y^ - Y^- + Y^- G dn| r < oo)p (3.18) 
= H*G{du), (3.19) 

where G{du) is the conditional distribution P(y.^- G du\ r < oo). From corollary 4.6 of Huzak et 
al. (2004), G{du) is equal to unconditionnal distribution P(yT— G du). Corollary 4.10 of Huzak 
et al. (2004) implies that 

P(y^- G du) = P(sup(-ct - aBt) G du), (3.20) 
t>o 

from which it immediately follows that G = G. Thus we have proved the the statement 2 of 
proposition. ■ 

For the right inverse of we shall write (f) on [0, oo). That is to say, for each > 0, 

$(g) = sup{a > : V'(a) = q}- (3.21) 
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Note that the properties of ■0 for such a Levy process X, imply that (f>{q) > for q > 0. Further 
0(0) = 0, since V'(0) = c - E[Si] > 0. 

Let us define the probabihty measure P by the density process 

= e-*('?)^*-9* (3.22) 

where ^{q) is the right inverse of ip defined in (|3.2ip . Note that under P the process Y introduced 
in (|2.ip . is still a spectrally positive Levy process, and still drifts to —oo [see Kyprianou (2006)]. 
In addition, the process Y keeps the same form under P and then, Yf = —ct + St — crBt, where 
c = c + (T^$(g) and Bt = Bt — a^{q)t are respectively the premium rate and the standard 
Brownian motion under the probability measure P. 

We denote by u the Levy measure of 5 under the change of measure P and then, 

?(dn) = e-*(«)"j^(d'u). (3.23) 

Let 

nS,] /(o,oo)^e-^('')V(d,) 

then by (12. 4p . < p < 1 and the net profit condition is well preserved under P. Note that 
the distribution function of — inft>o(.^t) = sup^>o(— ct — aBt) under the change of measure, is 
denoted by G. Recall from (j3.6p that the density G is given via its Laplace transform by 

c/3 



oo ^ 

e-^yG{dy) = — /3>0; (3.25) 



and then, G is given explicitly as an exponential distribution function with parameter 2c/ a"^. 
We denote by H the distribution of 1^ — Y^- under the probability measure P, i.e 

H{du) = ]P(y^ - F^- G ; r < oo), for li > 0. (3.26) 

Moroever, since the characteristics of the risk process are preserved under the change of measure, 
we can derive a result analogous to Proposition 13.11 under P as: 

Proposition 3.2 Let Y he a spectrally positive Levy process defined in i2.3\) . 

1- The distribution ofYj- — Y^- on the set r < oo under the probability measure P is given by 

1 ^ 

H{du) = - v{du + y)dy 
c Jo 

„-'S>{q)u i-oo 

= / e-*(^)^z/(dn + y)dy ; u > 0. (3.27) 

c + <^{q)a^ Jo 

2- The distribution of Yj- under the probability measure ¥ is given by 

F{Y^ edu; T <oo) = H* G{du) ; li > 0. (3.28) 

In the following paragraph, we shall give the ruin probability under the change of measure 
defined by ([322]). _ 

We denote by 6[x) the ruin probability under the probability measure P, that is 

e{x) = P(sup Yt<x) , X > 0. 

For g > 0, the following proposition gives the Pollazek-Hinchin formula for the survival proba- 
bility under the change measure P. 
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Proposition 3.3 The survival probability of the general perturbed risk process introduced in 1^2. 2\l 
is given by 

oo 

l-e{x) = {l-p)Y, (^*^"^ * ; (3.29) 

n=0 

where ^ 

L{dy) = ^z7(y,oo)dyl{y>o} , 

and G{-) is an exponential distribution function with parameter Icja^ , /*" (n > 1) denotes the 
n-fold convolution of f with itself and f*^ is the distribution function corresponding to the Dirac 
measure at zero. 

Proof. Using a similar method to that in Huzak et al. (2004), we have by taking hmits 
in the Laplace transform of the infimum evaluated at an independant exponential time Cq with 
parameter g' > [see chapter VIII in Kyprianou (2006)] , 

= ^'{Q+)J- for /? > 0. (3.30) 

V'(/5) 

Let us now compute (j3.30p in terms of p, L and G, where p\s the parameter given by (j3.24p . 
and L and G are respectively the Laplace transforms of L and G. Equation (j3.30p is equal to 

i ^^^K = {l-p)hp)f^[pG{P)L{P)T- (3.31) 

^ 1 - pG{^)L{^) n=0 
By inverting the Laplace transform (I3.3ip . we obtain 

oo 

P(Foo <x) = l-eix) = {l-p) ^ p"(G*("+i) * Z™) (x). (3.32) 

n=0 

■ 

We now introduce the so-called q-scale function {W^'^\ Q > 0} of the process X. For every 
q > 0, there exists a function VK^^) • — s, [0,oo) such that W^'^\y) = for all y < and 
otherwise absolutely continuous on (0, oo) satisfying 

/•CO 1 

/ e-^yW^'^\y)dy = —- , for X > ^q) , (3.33) 

Jo V'(A) - q 

where ^{q) is the largest solution of 'ip{/3) = q defined in (|3.2ip . 

For short, we shah write W^°^ = W. Let us introduce the 0-scale function under P, which 
we write as related to the q^-scale function of X, W'-i^ via the relation 

iy('?)(y) = e*('')^VF^(,)(2/). (3.34) 

The reader is otherwise referred to Bertoin (1996) and Chapter 8 of Kyprianou (2006) for a fuller 
account. 

In the following, we shall describe the discounted penalty function introduced in Definition 
12.21 in terms of densities G, H and g-scale function of the spectrally negative process X which 
we have briefly introduced above. 
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4 Extension of the Expected Discounted Penalty Function 

At this point, we recall that the main objective of this paper is to write an expression for the 
extended EDPF in (I2.1ip . But before we can write out such an expression we need one more 
intermediate result that has to do with the change of measure defined through the density process 
in (I322D- 

Let us start by giving the standard EDPF defined in (|2.6p in terms of the convolution product 
of two functions depending on Levy measure v and g-scale function which is defined by (13.33^ . 
Recall that u; is a non-negative bounded measurable function on M x M such that 0) = 0. 

Lemma 4.1 Consider the risk model in HI 



1. The distribution of the overshoot at Tx, Tx, is given by 

Tx{du) = / / e-*(5)V(dn -x + v)WLAx - y)dydv , (4.1) 
Jo Jo 

where Tx{du) = P(i^^ S du] Tx < oo). 

2. For q>0, the EDPF (j){w,q,x), as introduced in Definition \2.1[ is given by, 

h * f2{x) , (4.2) 



where 



W^^^x) - <^{q)W^''\x) , (4.3) 



and 



poo p 

f^{x) = e*(«)M e-*(«)M w{u,v)u{du + v)dv , (4.4) 

Jx J{0,oo) 

for u > x,v > and < y < x Av. 

Proof. 1). Follow the method used in BifRs and Kyprianou (2010) [see also the end of 
Section 8.4 of Kyprianou (2006)] and recalling that X drifts to oo (and hence $(0) = 0), we 
know that 

F{Yr^ - X e du,x - y.- G dv; Tx < oo) = z/((in + v)[W{x) - W{x - v)]dv 

pv 

= v{du + v) I W'{x — y)dydv , 

(4.5) 







for > 0, f > and < y < x A f , where W' is a version of the density of W . 
Using the previous equality under the change of measure P, we obtain the identity 

'^{Yr^ — X G du, X — y_- G dv ; Tx < oo) = u{du + v) 1^4(g)(^ ~ y)dydv 

= e-*('?)("+^)i/((iu + ^;) 

pv 

X / W^{g)ix - y)dydv , 



(4.6) 
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for u>0,v>0,0<y<xAv and then 



E 



e '^'^'^wiYr^ - x,x -Y - );Tx < oo 



E 



„'I>(g)(w+x) 



(0,00) J (0,00) 
xF^Yr^ — X e du,x — y.- G dv; Tx < 00) 



J(0,oo) 



where 

Mx) = e*(«)^- 
By using (13.34^ we get, 



h * f2{x) , 



/ w{u,v)i^{du + v)dv. 

J{0,oo) 



(4.7) 
(4.8) 



(4.9) 



from which, statement in 1) holds. 

2). From (j4.6p . we deduce (j4.ip by writing, 

Tx{du) = "Fi^Tx — X G du — x; Tx < 00) 
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for u > x,v > and < y < x A f . I 



"*(«)"z^((in -x + w)VF4(^)(x - y)(i?/d'(; 



(4.10) 



Remark 4.1 Equation J^.^p is equivalent to the following equality which describes the Gerber- 
Shiu function in terms of scale function of risk process; 



) = ! [ e-*(^)(^-^)u;(n,?;)i/(dtx + t;)[t^*(«)(x) -Ty*(5)(x-7;)]d7; 
J {0,00) J {0,00) 

00 roc 

w{v,u)[e-^^i>W^'i\x) - W^i\x - v)]i^{du + v)dv , (4.11) 



[w,q,x 



Jo 

by using ^3.34^ , where w is a bounded measurable function such that w{-,0) = 0. 

Identity (|4.1ip is given in more general form in Biffis and Kyprianou (2010), where the EDPF 
also includes the size of the last minimum before ruin x — Y-r^-. 

Recall that from Section [3] that N is the random number given by (|2.1Up . We give in the 
following proposition the distribution of A'^ under the probability measure defined by ()3.22p . 



Proposition 4.1 The distribution of N on {tx < 00} is given by 

F{N = n, < 00) = (1 - p)r{l - 0{x)) ■ 

where n > 0. 



(4.12) 
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Proof. By the strong Markov property of Y under P, we can identify the distribution of N on 
{tx < 00} as 

P(iV = 0, Tx < 00) = P(r(^) =oo\tx< oo)P(r^ < 00) 
= P(r = oo)(l-^(x)) 

= (4.13) 

where in the last equality we have used Corollary 4.5 of Huzak et al. (2004) under the probability 
measure P. 
For n > 1, 

¥{N = n,Tx< 00) = P(r("+i) = oo|r(") < oo)P(r(") < oo|t(""^) < 00) x 

... X ]P(t(^) < oolr^; < (X))P(t^ < 00) 
= P(r = oo)P(t < oo)"(l - e{x)) 

= (l-p)p"(l-^(x)). (4.14) 



Now, we are finally in a position to present the first and main contribution of this paper. 
In the following theorem, we give an expression for the extended EDPF P{F; q; x) given in 
Definition 12.21 This expression is given in terms of the g-scale function. Levy measure, and the 
densities G and H which are introduced in Section [3j Recall that F = (-F„)n>o is a sequence of 
non-negative measurable functions from M+ x ]R+ to M such that Fq{.,x) = 0. 

Theorem 4.1 Consider the risk model in \2. S]) . For q >0, the EDPF P{F, q, x), as introduced 
in Definition \2.S\ is given by, 



00 ^ ^ 

/ e*(«)("+-)F„+i(t;,n + ^;) 

1 J(x,oo) J (0.00) 



[w,q,x 

n=l 



' (a;,oo) J (0,00) 

H * G{du)H*'' * G*" * fx[dv) , (4.15) 



where w is a mesurahle Borel- function satisfying w{u,v) = Fi{x — v,u — x) for u, v >0 and 



(w, q,x) = E 



e~''^^ w{x - y - ,Yr^ - x) ;Tx < 00 . (4.16) 



Proof. We prove the result in three steps. Let us suppose F = {Fn)n>i is a sequence of 
non-negative measurable functions from x to M, x and q > 0. 
Step 1: We prove in this step that 

P(y,(fc) - y,(.-i) e d2/|T« < 00) = Ih* G{dy)- (4.17) 

where k > 2 and y > 0. The proof of ()4.17p follows by an application of Theorem 4.7, Huzak et 
al (2004). Recall that under the measure change P, Y has the same form and it can be written 
as Yt = —ct + St — cBt- By using the Markov proprety of Y at t^^'^ 

P(y^(fc) - y^(fe-i) G dy|rW < 00) = P(y, G dy\T < 00) 

= LH*G{dy), (4.18) 
P 
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where in the last equahty we have used Equation (I3.28P and the identity P(r < oo) 
Corollary 4.6, Huzak et al. (2004)]. 
Step 2: Next we prove that 

P(y,(.) € dy|r(^) < oo) = \ (H*' * G*^ * f,){dy), 

p'^yl — v[x)) 

where k >2 and y >0. 



p [see 



(4.19) 



P(y,w G dylrW < oo) = P(y,(.) - Y^,,-,, + ... + - Yr^ + Yr^ G dy|T^^^ < oo 



(fc) 



1 (i/*W*G*W*r,)(dy) \ 



1 - 6'(x) 



(4.20) 



by using the independent increments of Y and ()4.17p . where Tx{dy) = P(1V^ ^ dy ]Tx < oo). 

Step 3: Before using the conclusions of step 1 and step 2, let us write P{F,q,x) as an 
expansion under P, 



oo n 



P{F,q,x) = J^E[j^e-5-''Vfc+i(y^w,y,(.+i));iV = ?i;r, <oo 

n=0 fc=0 

= E [e-"^- Fi (y - , ; iV = 0; < oo) 

oo n 



+ 2^ E ^ 2^ e"' ^^^+1 (^r{^) , ^C^+i) ) ; iV = n; < oo 

n=l fe=0 



n=l fe=0 

E [e-"^- Fi (y - , y^J ; iV = 0; < 

oo 

+ (e [e-^^- Fi (y - , y^J ; iV = n; < oo 

n=l 
n 

+ ^ E [e-«^''^'^ Fk+i (y,(fe) , y,(..+i) ) ; TV = n , < oo 



k=l 



where in the last equality we have used Proposition 14.11 Then P{F,q,x) is equal to 



T:r: < OO 



(i-m-oix)) 



E[e*(«)^-Fi(y-,y,^ 

oo n 

+ ^ (E[e-''-^Fi(y-,y.j|T, < oo)] +^E[Ffc+i(y,(.),y, 



(4.21) 



{.+!)) Ir^'^+i) < oo 



1.22) 
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oo 

E[Fi(y y.j|r, < oo)] (1 - e{x)){i - + ^p-] 

n=l 

oo n 

+ Y.(Y.^[e'''^'^^ri^+"'' Fk+, , ) I T^'^+i) < ^] ) ( 1 - ( 1 - ^(x) ) 

n=l fc=l 

oo n „ „ 

'-^-'(0,00) J{0,oo) 



P(F^(..+i) - y^(fe) G dn.y^tfe) € dz;|T('=+^) < oo) (1 - p)p"(l - 0(x)) 



(4.23) 
(4.24) 



1 JfCoo) JfO.oo) 



xP(y^(„+i) - y^(„) e r("+i) < oo)P(y,(„) e (i^;|r("+i) < oo))p"+Hl - ^(x)) 



(4.25) 



w,x) + y^ / e*(«)("+^)i^„+i(7;,^z + 7;)(^*G)(d^z) 

„=l^(0,oo) y(0,cx>) 



*r^)((i't;) 



where in the last equahty we have used ()4.17p and (j4.19p . 



1.26) 



Remark 4.2 Let F = {Fn)n>i he a sequence of hounded measurable functions. 

If F satisfies Fi{u,v) = w{x — u,v — x) and F„ = for n > 2, where w is a measurahle 
function such that w{-,0) = 0, then P{F,q,x) reduces to the classical EDPF defined by ^4-^ (^'^d 
then. 



P{F,q,x) 



[w,q,x) 



K 



"'?^-u;(x - y - , - x); T^. < oo 



t.27) 



which is completely characterized in terms of Levy measure and scale function in Lemma \4.1\ and 
Remark \4-l\ 



5 Capital injections 

In this subsection, we introduce the Expected Discounted Value of Capital Injections (EDVCI), 
which are necessary to keep the reserve process R above 0. In our context, if R goes under 
by jumping, we must apply control to prevent the process staying in ( — oo,0). In fact, we 
should inject capital only when the risk process becomes negative and only when, the new 
record infimum under (undershoot) is reached by a jump of a subordinator. 
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Recall that in Theorem 14. ll we have identified the extended EDPF P{F; q; x) defined by (12. 2p 
in terms of the g-scale function, Levy measure, and the densities G and H introduced in Section 
m Using the connection with the result of Theorem 14. H we will characterize the EDVCI for 
the risk process defined in (j2.2p . We will study more explicitly the classical case driven by the 
Cramer-Lundberg risk model [see Einsenberg and Schmidli (2011)]. 

5.1 Expected discounted value of capital injections (EDVCI) 

Recall that R is the risk process defined in (|2.2p by 

Rt:=x-Yt, t>0, (5.1) 

where Yt = —ct + St — crBf. 

We denote by Ct the cumulative capital injections up to time t. The controlled risk process 
is given by 

Rf ■.= x-Yt + Ct , t>0. (5.2) 
We have to inject the first capital when the risk process falls below zero. Let 

t(i) =t^ = inf {s >0,Rs < 0} 

denote the time of first ruin, and 



Ci — Rj 



X 



denote the first injection. At time r*-^-* the controlled risk process R'~" starts with initial capital 
equal to zero. 
For n > 1, let 

t("+^) =inf{t >r("),^,- <Rt} 
denote the time of the n-th injection. The size of the n-th injection becomes 

Cn = Rr^n) — R^(n+l) 

= Y^{n+i) — Y^(n), for n > 1. 

The accumulated injections can be described as 

N 

Ct = ^ C„l|^(„)<^|. 



(5.3) 



Let us define the expected discounted value of capital injections (EDVCI) as 



TV 

V{q,x)=E[Y,e-''-'"'''Cr, 



n=0 



Let introduce K{q, x) as 



and 



K{q, x) = E 



e '?^^;Ta; < oo 



ipiq,x) = E e 5^-(y^, - x);r^ < oo 



(5.4) 



(5.5) 



(5.6) 
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From Lemma 14.11 we can give an expression for (15. 5p and (15. 6p in terms of g-scale function and 
then 



and 



where 



ip{q,x) = h*h{x) (5.7) 



i^{q,x) = fi*t{x); (5.5 



h{x) = e*(^)^ / e-*(«)'' / uuidu + v)dv 

Jx J{0,oo) 



r>oo 

I I _ v)i^{du)dv, (5.9) 

' (i',oo) 



/■oo 

t(x) = e*('')M e-*('')V(t;,cx))(it;. (5.10) 



and 



Recah that in Theorem 14.11 we have identified the extended EDPF introduced in Definition 
2.21 in terms of the g-scale function, Levy measure and notions introduced in Section [3] for 
the model (j2.2p . By using connection with Theorem 14.11 the foUowing theorem is an exphcit 
characterization of the EDVCI defined by (15. 4p . This would be the second main contribution of 
this paper that extends similar results in Einsenberg and Schmidli (2011). 



Theorem 5.1 The Expected Discounted Value of Capital Injections EDVCI introduced in 1^5. 4\ ) 
is given by 

Viq, x) = ^{q, x) + ^^If^ k{x, a) , (5.11) 
where C(^,o") and 5{q,a) are given, respectively, by 

= (l + _lM?i^)[l (5A2) 

and 

2c r 2g i 

" ^iq){2c + ^{q)a^) U(g)(2c + ^{q)a^) + ^ " ' ^^'^^^ 

Equation (j5.1ip gives an explicit formula of the expected value which should be injected at each 
deficit time at and after ruin that will allow the insurance company to survive and continue its 
operations when the risk process continuous to jump downwards. 

Proof. Let us consider the sequence of functions F = {Fn)n>i as defined above. Since 
we suppose Fi{v,u) = u — x and Fn{v,u) = u — v for u > 0, v G M and n > 2, then the 
extended EDPF associated with F and q, P{F,q,x), defined in (|2.2p . is equal to V{q,x). By 
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using Theorem 14.11 we can easly derive (15. IIP and then 

+ 



^/ / e*(«)("+'')'u(^*G)(dn)(^*"*G*"*ri''))(d?;) 

n=l (^1°°) (0,oo) 

= 'f{q,x)+ [ e'^^i^^'u H * G{du) 

J{0,oo) 

oo „ 

= ip{q,x)+ e'^^'^^^'u H * G{du) e^^'^^''f^{dv) 

J(0,oo) 

^ oo ^ 

= ip{q,x)+ e'^^'^^^'u H * G{du) e'^^'^^f^idv) 



i(0,oo) 



(5.14) 



Recah that P{Yt-^ ^ dv Tx < oo) = Tx{dv), for > x, hence by (|3.22p . (II) is equal to 

E[e'^ii)y^.- Tx<oo] = E[e-«^^ < oo] 

= K{q,x). (5.15) 

(III) is equal to 

[ e*(«)"^*G((iu)= /" e*(«)"#(dn) / e^^")" G(dn) , (5.16) 

where 



f oo 

z.(y, oo)(iy. (5.17) 



c + $(g)(T2 7o 
Recall that 

V'(^(q)) = q , (5.18) 

by using integration by part, (|5.18|) is equivalent to 

C^iq) + ^ - ^(g) / e-*^'?)^ i/(y, oo)(i?/ = q (5.19) 

2 Jo 



and then, 



/■oo 

/ e-*(«)^z^(y,oo)dy = c + 



(5.20) 
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We have 



L 



(0,oo) 



G{du) 



2c 



re o-^ 



dy 



(o,oo; 

1 + ^ ...NO = a(g, 0-) 



2c + (t){q)a'^ 

By substituting (|5.20p in (j5.16p . we conclude that 

^(0,oo) 



(5.21) 



a(g,cr)[l 



$(g)(c+4$(g))- 



(5.22) 



(I) is equal to 

e*(9)«n H * G{du) 

(0,oo) 



e*(9)«[/ (n + r;)e*(«>G(d7;)] iJ(d- 



(0,oo) 



(0,oo) 



ei>('?)«[/ (n + r;)e*(«>G((i7;)] F(d- 



(0,oo) 



(0,oo) 



^<i>{q)u 



(0,oo) 



[n/ 



2c „f2c 



(5„$(<j))„ 



+ 







°° 2c _(2c_^( 



(it;] i7((it 



e*(^)"[a(g,a)(u + 



(0,oo) 

a(g,cj; 
+ 



a 



2c + ^>(g)f72 

2 



)] 



c + «>(g)a2 L2c + «>(g)a2 j(o „^) 

ue 



J(0,oo) 
1 



a(g,<T)- 



0" 



C + ^'(g)0-2 L2c + ^>(g)0-2 



(0,oo) 



''^^'^^y v{y,oo)dy 



v{y,oo)dy 



a(g, a) 



1 



1 



c+^>(q')fj2 L^2c + ^>(g)cr2 $(g) 
e-*('')^ ^(y,oo)dy+-i^E[5i] 

(0,oo) ^('7) 



) 



(5.23) 



Using (|5.20p . the last equality can be written as 



-2c 



c + $(g)a2 L$(g)(2c + ^>(g)CT2) 



c + 



a2$(g) g 



+ 



(/>(g)^ ^(g) 



2q 



^{q)ic + ^{q)a^) L$(g)(2c + $(g)a2) 
5(g,cr) 



+ P-1 



(5.24) 
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and then (I) is equal to 

[ e'^^'i^''uH*G{du) = 5[q,a). (5.25) 

i(0,oo) 

In addition, by using the identification above of (II) and (III), (|5.14p is equal to 

oo 

(^(g,x) + (5(g,cr) ^K{x,a)i{q,aY 

n=0 

= (^(g,x) + -^%^K(x,a) (5.26) 

and the theorem follows. ■ 

In the next subsection, we illustrate the previous result by a specific example of Cramer- 
Lundenberg risk model when the Brownian component vanishes. 

5.2 Classical risk model 

Let us consider the classical risk model. The number of claims is assumed to follow a Poisson 
process {Nt)t>Q with intensity A. We denote by {Zn)n>i the claim sizes which are independant 
of {Nt)t>Oi positive and iid with distribution function K and first moment /u. The aggregate 
claim process is given by 

Nt 

Rt = x + ct-^Zi, (5.27) 

i=l 

where x is the initial capital and c is the premium rate. 

The surplus process introduced by (15. 2p then has the form 

Nt 

Rf := x + ct-^Zi + Ct , t>0 . (5.28) 

i=l 

Note that St = ^^^i Zi and then E[5i] = Xfi. However, the net profit condition (12. 4p becomes 
c > Xn and then p = A/x/c. Since cr = 0, it follows from Theorem 15.11 that the EDVCI for 
classical model introduced above is given by 

Viq,x) = ^{q,x) + j^^^^t,{x,0), (5.29) 



where 



1 . n^ 1 



^(^'0) = 1 - ^TtV' and%,0) = --f- --f- + p-l =- (5.30) 

$(g)c ^{q) l^{q)c J c$(g)2 



(c - XfJ.)<^{q) 



Since the Brownian component vanishes, y is a process with bounded variation and then, from 
Hg){dx) = T.n=0'- 



Kyprianou (2006), W^,g){dx) = T.n=oV*"'{dx), where 



r](dx) = -u(x,oo)dx 
c 



e-*(9)«i,((iu) dx. (5.31) 

{x,od) 
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Note that Equation (14.90 reduces to 



. oo „ 

/i(^) = -E / (5.32) 

when in the last equahty, we have used the identity u{du) = \K[du). In addition Equations 
dSj]) and (|5lT)l) reduce to 

h{x) = -e^ii> r e-^^'i^'' I {u-v)K{du)dv (5.33) 

and 

t{x) = -e*(«)^ / e-*('?>(l - F{v))dv. (5.34) 



Consequently, Equation 15.291 reduces to the expression of the EDVCI for classical model (|5.27p 
given in Einsenberg and Schmidli (2011). 



6 Conclusion 

In this paper we have generalized the Expected Discounted Penality Function (EDPF) indro- 
duced by Gerber and Shiu (1997, 1998) to include the successive minima reached by the risk 
process because claim after ruin. In addition to the surplus before ruin and the deficit at ruin, 
we have added to the EDPF the expectation of a sequence of discounted functions of minima in 
the context of subordinator risk model perturbed by a Brownian motion. By using some results 
in Huzak et al. (2004) and developments in theory of fluctuations for spectrally negative Levy 
processes, we have derived an exploit expression of this extended EDPF. 

Our generalization of the EDPF includes information on the path behavior of the risk process 
not only in a neighborhood of the ruin time, but also after ruin. Such information is relevant for 
risk management process aimed at preventing successive occurences of the insolvency events. In 
addition to the classical EDPF introduced by Gerber and Shiu (1997, 1998), the new extension 
of EDPF contains a sequence of expected discounted functions of successive minima reached 
by jumps after ruin. This sequence of EDPF has many interesting potential applications. For 
example, it could be used as a predictive tool for successive deficit times after ruin. In particular, 
we have used this extended EDPF to derive explicitly the Expected Discounted Value of Capital 
Injections EDVCI which are necessary to keep the risk process above zero. 

Inspired by results of Huzak et al. (2004) and developpements in fluctuation theory for 
spectrally negative Levy processes, we provide a characterization for this extended EDPF in 
a setting involving a cumulative claims modelled by a subordinator, and spectrally negative 
perturbation. We illustrate how the ESDPF can be used to compute the expected discounted 
value of capital injections (EDVCI) for Brownian perturbed risk model. The main contributions 
of this paper are found in Theorems 14.11 and 15.11 These two results give expressions for the 
extended EDPF given by Definition 12.21 and the EDVCI introduced in Subsection 15.11 These 
expressions can be easily computed by considering particular examples of subordinators. 

In oder to make these results as explicit as possible, we have used the risk model driven by 
Brownian perturbed subordinator process. Further work is needed in order to give these results 
for more general risk processes driven by a spectrally negative Levy process. 



20 



References 

[I] Asmussen, S. (2003). Applied Probability and Queues, 2nd edition. Springer- Ver lag, Berlin, Heidel- 
berg, New York. 

[2] Asmussen, S. and Albrecher, H. (2010). Ruin Probabilities. World Scientific Publishing, London. 
[3] Bertoin, J. (1996). Levy processes. Conibridge university press. 

[4] Ben Salali, Z. and Morales, M. (2012). Levy Systems and the Time Value of Ruin for Markov Additive 
Processes . European Actuarial Journal, 2(2), 289-317. 

[5] Biffis,E. and Kyprianou, A. (2010). A note on scale function and the time value of ruin. Insurrance: 
Mathematics and Economics. 46, 85-91. 

[6] Biffis,E. and Morales, M. (2010). On a generalization of the Gerber-Shiu function to path- dependent 
penalities. Insurance: Mathematics and Economics. 46, 92-97. 

[7] Dufresne, F. and Gerber, H. (1991). Risk theory for compounted Poisson process that is perturbed by 
diffusion. Insurance Math.Econom.lO, 51-59. 

[8] Dufresne, F., Gerber and Shin, E. W. (1991). Risk theory with gamma process. Astin Bull. 21 177-192. 

[9] Doney, R. and Kyprianou, A. (2006). Overshoots and undershoots of Levy processes. Ann. Appl. 
Probab. 16(1), 91-106. 

[10] Einsenberg, I. and Schmidli, H. (2011). Minimising expected discounted capital injections by rein- 
surance in a classical risk model. Scand. Actuarial J. 3, 155-176. 

[II] Huzak, M., Perman, M., Sikic, H. and Vondracek, Z.. (2004). Ruin probabilities and decompositions 
for general perturbed risk processes. The Annals of Applied probability. 14(3), 1378-1397. 

]12] Furrer, H. (1998). Risk processes perturbed by a-stable Levy motion.. Scand. Actuar. J. 1, 59-74. 

]13] Garrido, J. and Morales, M. (2006). On the expected penality function for Levy risk process.. N. Am. 
Actuar. J. 10(4), 196-218. 

]14] Gerber, H. and Shiu, E. (1997). The joint distribution of the time of ruin, the surplus immediately 
before ruin, and the deficit at ruin. Insur. Math. Econ. 21, 129-137. 

]15] Gerber, H. and Shiu, E. (1998). On the time value of ruin. North Am. Actuar. J. 2(1), 48-78. 

[16] Gerber. H. and Shiu, E. (1998b). Pricing Perpetual Options for Jump Processes. North Am. Actuar. 
J 2(3), 101-112. 

]17] Gerber. H. and Landry. B. (1998). On a discounted Penalty at Ruin in a jump- diffusion and Perpetual 
Put option. Insur. Math. Econ. 22, 263-276. 

]18] Kyprianou, A. (2006). Introductory Lectures on Levy Processes with Aplications . Universitext, 
Springer. 



21 



